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Recently, a framework for analyzing time series by constructing an associated complex network has attracted 
significant research interest. One of the advantages of the complex network method for studying time series is 
that complex network theory provides a tool to describe either important nodes, or structures that exist in the 
networks, at different topological scale. This can then provide distinct information for time series of different 
dynamical systems. In this paper, we systematically investigate the recurrence-based phase space network of 
order k that has previously been used to specify different types of dynamics in terms of the motif ranking 
from a different perspective. Globally, we find that the network size scales with different scale exponents and 
the degree distribution follows a quasi-symmetric bell shape around the value of 2k with different values of 
degree variance from periodic to chaotic Rossler systems. Local network properties such as the vertex degree, 
the clustering coefficients and betweenness centrality are found to be sensitive to the local stability of the 
orbits and hence contain complementary information. 



Methods of nonlinear time series analysis seek to 
quantify deterministic structure in a dynamical 
system from measured time series data. To date, 
these methods are based on dynamical systems 
theory and the indices of interest quantify fea- 
tures related to the (possible) underlying deter- 
ministic dynamics. Measures such as correlation 
dimension and Lyapunov exponents are widely 
applied and have found utility in a diverse range 
of sciences. The underlying principle of these 
methods is that the various embedding theorems 
provide conditions under which scalar time se- 
ries measurements can be transformed into a ge- 
ometric object which (under suitably generic con- 
ditions) is equivalent to the (presumed) original 
attractor. Recently, a new set of transforma- 
tions have gained considerable interests. With 
these transformations it is possible to represent a 
scalar time series as a complex network. Previous 
work has already shown that these complex net- 
works retain certain features of the original dy- 
namical systems and therefore the arsenal of com- 
putational complex network science and theoret- 
ical graph theory may now be applied to analyse 
structural properties of time series data — and by 
implication the underlying dynamical system. In 
this paper we provide the first thorough analysis 
of quantitative measures of the complex network 
properties of the recurrence-based phase space 



networks. We study a range of network based 
measures at the local-, meso- and macro-scale and 
show how these measures relate to the dynamical 
properties of the underling system. We also ex- 
tend previous work by providing a rigorous com- 
putational study of the effect of noise in the time 
series data on the results: both for simulated and 
experimental data. 
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I. INTRODUCTION 

Time series analysis and its applications have long 
been important for the understanding of complex sys- 
tems in many research fields such as sociology, biology 
and society^. Early studies mainly focused on the time 
series and its phase space embedding and this led to non- 
trivial features estimated from data: long-range correla- 
tions, scale invariance and so on. Meanwhile, a lot of 
recent work has shown that natural and artificial sys- 
tems in diverse fields can be also described with complex 
networks^Sl which are composed of a number of nodes 
that are interconnected by a set of edges. Although 
time series and complex networks share many relevant 
applications, the bridge between time series and complex 
network has not appeared until Zhang and SmalP pro- 
posed the method of transforming pseudoperiodic time 
series into the so-called cycle networks in 2006. Since 
then, many methods for transforming time series into 
complex networks have been designe d, suc h as recur- 
rence networks^HH, correlation network d 10 * 11 ! and visibil- 
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ity graphs^. By going from time series to complex net- 
works, it is possible to understand the correlation struc- 
ture and dynamical properties of the time series by using 
the mathematical knowledge from graph theory and com- 
plex network science. 

In analyzing complex networks, one of the greatest ad- 
vantages is the existence of multiple scales of identifica- 
tion and description of the important nodes or structures 
that exist in the networks. 

For a given a network, even simple global statistics 
can provide a profound insight into the network itself. 
The most intensively studied "macroscopic" concepts are 
the average path length, average clustering coefficient, 
and degree distribution. Examples include the two well- 
known characterization of scale-free and small world fea- 
tures widely found in networks. The former was based 
on the observation that the degree distributions of many 
real networks have a power-law form. The latter refers to 
the small world phenomenon which is found especially in 
many social networks and is characterized by a low aver- 
age path length and high average clustering coefficient. 

When investigating complex networks at the "micro- 
scopic" scale, local properties of the network can also 
be important. Identifying the most vital nodes has also 
attracted academic research interest since these nodes 
could play a crucial role in the organization of the net- 
work and thus has significant impact on the dynam- 
ics of networks, such as their synchronization, epidemic 
spreading and so on. Th e prop osed network statistics 
to measure node centrality'^ include degree central- 
ity, local clustering coefficient, betweenness centrality, lo- 
cal subgraph centrality, eigenvector centrality, closeness 
centrality, and so on. Though considered as local proper- 
ties, these network statistics actually capture the network 
properties on a different scale. For example, the degree 
centrality which is defined as the number of neighbors of a 
vertex, measures only the local connections for the given 
node. However, the local clustering coefficient reflects the 
importance of the node in the network at a more macro- 
scopic level than degree centrality since it measures the 
connectivity between the neighbors of a node. Similarly, 
the betweenness centrality of a node involves calculating 
the shortest paths between pairs of vertices that pass the 
node on a graph and therefore measures the importance 
of a node at an even more macroscopic level. 

Moreover, much recent research attention has been put 
into studying the structural important features of the 
network, such as the motif, cyclic, assortative patterns, 
and other community structures"22l. These "meso- 
scopic" structures which consist of both nodes and edges 
play an important structural and functional role in com- 
plex networks. 

The basic idea of this paper is that the conversion from 
time series to complex networks provides a multiscale 
characterization originating from the network analysis, 
which in-turn gives a deeper insight into the dynamical 
origin of the time series. The key questions are what 
information of the time series is contained inside the net- 



work representation and how may such multiple scales of 
identification and description be related to the time se- 
ries. To address these questions, much research effort has 
been devoted to networks based on the important concept 
of recurrences from time series in the phase space. Tradi- 
tionally, to construct a recurrence network following the 
idea of a recurrence plot proposed by Eckmann et al™ — 
the time series is mapped into a set of points in the em- 
bedded phase space, each point represents a node of the 
network, and two nodes are linked together when their 
phase space distance is smaller than a selected thresh- 
old e. By choosing an appropriate e, the local phase 
space properties can be best preserved^. These so-called 
epsilon recurrence networks are phase space density de- 
pendent, and link the local properties of the network 
with the phase space properties of the attractor directly. 
This approach allows one to uncover similarities between 
the phase space structures of various complex systems, 
which can provide deeper insights into the complex sys- 
tems. In 2008, Xu et al. suggested the so-called phase 
space networks constructed by linking each node with its 
k nearest neighbors^ in the reconstructed phase space. 
By studying the occurrence frequencies of the connected 
motif patterns, this methods can be used to distinguish 
time series with different dynamics by their so-called su- 
perfamily phenomenon. An advantage of this approach 
is that the phase space networks are topologically invari- 
ant, since they depict the relative proximity relationships 
among a fixed number of neighboring points in the phase 
space attractor. 

As in Xu's phase space networks (Ref. the net- 
work motifs which are the building blocks of the complex 
networks reflect different local structure properties' at 
the microscopic level. Besides the motifs, however, other 
characteristics which are associated with local and global 
properties of the network and characterize the impor- 
tance of a node in the network from different scales have 
been left unstudied. In this paper we address this defi- 
ciency. To do so, we will focus exclusively on the phase 
space network of Ref. [24] and not the recurrence meth- 
ods of Refs . 171 lr?l and [231 For an extensive comparison of 
these methods, and the other emerging alternatives, see 
Ref. [Ml 
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FIG. 1. All six motifs of size four denoted by A, B, C, 
D, E and F. 

The reminder of the paper is organised as follows: 
we first quantify these recurrence-based phase space 
networ ks (u sing the adaptive nearest neighbor network 
methocP^H) v i a the standard topological statistics, in 
addition to the motif ranking, to show how such global 
and local network statistics can provide new informa- 
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tion about the phase space geometry beyond time series. 
Then we apply the same analysis to noisy experimental 
data from continuous musical tones produced on a stan- 
dard clarinet. In the last section of this paper we provide 
a brief conclusion. 



II. QUANTITATIVE ASSESSMENT OF 
RECURRENCE-BASED PHASE SPACE NETWORKS 

The recurrence-based phase space network representa- 
tive for a given time series is generated following the same 
method used in Ref . [24] and summarized in Ref . [26j The 
first step is to embed the time series in an appropriate 
phase space. After proper embedding, each point in the 
phase space is taken as a node of the network and will be 
assigned with k new nearest neighbors. Eligible k neigh- 
bors should exclude points within the same orbit as well 
as those which are already neighbors. On the one hand 
this ensures that the neighbors are spatial neighbors in 
the phase space rather than the temporal neighbors, on 
the other hand this avoids multiple links between two 
nodes and provides a measure of phase space density. By 
constructing network in this way, we are able to enforce 
a threshold adaptively according to the local inter-point 
density of the attractor and thus obtain a network rep- 
resentative of time series which inherits structure of the 
distinct local phase space properties from different dy- 
namical systems. 

For a low-period periodic time series, the associated 
network is regular as the points are arranged in an or- 
derly manner around an orbit in the phase space. The 
network forms a large loop shape and its size increases 
proportionately to the number of points of the embed- 
ded data. If we add noise to the data, the dimension 
of the dynamics increases while the homogenous distri- 
bution of the points remain unchanged: we still get the 
loop structure, but the diameter of the network decreases. 
For chaotic time series, the reconstructed phase space be- 
comes heterogenous and may have some fractal proper- 
ties and thus thus the structure of the associated network 
is more heterogeneous. 

As a typical example for a continuous dynamical sys- 
tem, we study the Rossler system. We make use of the 
data from the x component of the Rossler system by 
solving x' = — (y + z), y 1 — x + ay, z' = b + (x — c)z 
with a = 0.1 and b = 0.1, using the fourth-fifth or- 
der Runge-Kutta method (ODE45 in MATLAB) with a 
fixed sampling rate of h = 0.1. By tuning the parameter 
c, low-period periodic, high-period periodic and chaotic 
data can be obtained. Figure [2] shows the bifurcation 
diagram, with parameter c G [1, 18] divided into 1000 in- 
tervals. The data is mapped into a space of dimension 
d e = 10 that is sufficiently large to study the topological 
structure with the time delay r which is chosen to be 
the first minimum of the mutual information. To con- 
struct a network from the time series, we take each point 
in the embedded space as a node and link it with its 4 



nearest neighbors which will lead to a network with the 
mean degree equal to 8. The choice of 4 neighbors fol- 
lows Ref.|21]to guarantee a robust results — however, as 
noted in that reference, this choice does not significantly 
affect the results. Nonetheless, the number of nodes of 
the network is equal to N — r • (d e — 1) where N is the 
time series length. 
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FIG. 2. Bifurcation diagram illustrating the dependence 
of x component of Rossler system on the parameter c. 
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FIG. 3. Dependence of the average path length L(N) 
on the length N of the time series for the Rossler sys- 
tem in different dynamical regimes. Best fitting mod- 
els are L(N) = 0.0857V ' 96 and L(N) = 0.059iV 100 for 
the low-order periodic Rossler system of period 2 and 
3 respectively. In comparison, best fitting models are 
L(N) = 0.48iV ' 67 and L{N) = 0.78iV ' 43 for the chaotic 
Rossler system with c = 9 and c = 18 respectively. 



Global network characteristics 2. Clustering coefficient 



1. Average path length 



To investigate the distinction between qualitatively dif- 
ferent types of dynamics in more detail, we first study the 
average path length of the phase space network as a func- 
tion of the time series length. Figure [3] shows how the 
average path length of the phase space network changes 
with respect to the increasing length of the time series. 
Note that for shorter time series, the average path length 
can be affected by the choice of the initial point of the 
flow data, so the desired results are taken as an aver- 
age over several networks from different segments of the 
same parameter. We find that the average path length 
L(N) scales with N with different scaling exponents for 
low-order periodic and chaotic Rossler systems. The av- 
erage path length L(N) increases linearly with time series 
lengths N for low-order periodic time series with periods 
2 and 3, but exponentially for the chaotic Rossler system. 
For high-order periodic cases, we find transient behav- 
ior when the size N is small because the corresponding 
time series can be too short to cover enough information 
with respect to the periodicity. Moreover, the slopes for 
the average path length L(N) are larger than the two 
chaotic Rossler cases but smaller than the low-order pe- 
riodic cases due to the interactions between two nearby 
orbits appearing in the phase space. 



The global clustering is closely related to motif rank- 
ing. Motif F, which denotes a fully-connected subgraph 
of order 4 as shown in Fig. [T] indicates strong mutual 
coupling between nodes, and hence we may expect that 
the clustering coefficient for the periodic Rossler will be 
larger than the chaotic Rossler, because motif F occurs 
more frequently in periodic phase space networks rather 
than networks from chaotic data of the same Rossler sys- 
tem as mentioned in Ref.l2"4l From Fig. [4] we can see that 
the clustering coefficients are rather stable with respect 
to an increasing N and converges to around 0.68 for the 
low-order periodic regimes with periods of 2 or 3 as the 
network size increases. Similar behavior can be found in 
the chaotic Rossler system with c = 18. The clustering 
coefficients are also stable but converge to a lower value 
around 0.58. In comparison, significant transient effect 
can be found in the remaining three cases. The curve 
for the chaotic Rossler system with c = 9 first rises to 
the same high level as in the low-order periodic cases for 
the N between 1000 and 2000. Then, it decreases and 
finally converges to a low level similar to the chaotic case. 
The curve for the Rossler system of period 6 goes up and 
down and finally approaches the same value as the low- 
order periodic cases. Although the value of the clustering 
coefficient for the Rossler system of period 8 is closer to 
that of the chaotic cases, we may infer that, as the size 
of the network increases, the curve will finally converge 
to the low-order periodic values. 
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FIG. 4. Dependence of the global clustering coefficient 
C(N) on the length N of the time series. 



We also look into other global properties including 
the clustering coefficient and the degree distribution for 
Rossler systems of different dynamical type (chaotic, 
noisy or periodic motion). 



3. Degree distribution 

Degree distribution itself measures the global network 
connectivity and can therefore capture the global net- 
work topological properties inherited from the time se- 
ries. Figure [5ja) shows the average degree distributions 
for networks of several realizations for different initial 
conditions associated with 10000 points of the Rossler 
system with different dynamics. We find that the aver- 
age degree distribution follows a quasi-symmetric bell- 
shape for networks from the Rossler system in different 
dynamical regimes. The degree at which p(k) reaches its 
maximum is 8, which is exactly the average degree of the 
network (i.e. mode and median are equal). The peak val- 
ues of p(k) for networks generated from periodic data are 
all above 0.3 and the variance of the degrees with bell- 
shaped distribution are rather small, indicating that the 
networks are homogenous. In contrast, for the chaotic 
Rossler system, p(k) is lower than 0.3 and the variance 
of the degrees are larger. 

We plot the variance of degrees in the phase space 
network as a function of the bifurcation parameter c in 
Fig. [5jb) and find a much clearer relationship between 
the bifurcation diagram and the variance of degrees in 
the network. We make use of 10 realizations from ran- 
dom initial points. For each of the 10 realizations, there 
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FIG. 5. (a) The degree distribution for the Rossler 
system of different dynamical types, (b) Red line: The 
calculated periods according to different values of param- 
eter c. Blue line: The variance of degrees with respect to 
the change of parameter c. 



are 10000 points for every step of c. In Fig. [5jb) the 
red line shows the period while the blue line displays the 
variance of degrees with respect to the gradual changes 
of parameter c. To get the periods as a function of the 
parameter c, we first detect the locations of maximum 
points at cycle scale. If the difference of any two max- 
imum points is larger than a given threshold, we group 
them into different periodic orbits. Finally the number 
of orbits is considered as the period of the data. Due to 
the limit of numerical precision, if the numerical num- 
ber of periods exceeds a certain large value, then that is 
taken to be an indicator of chaos (that is, we are forced to 
consider aperiodicity and periodicity with period longer 
than some arbitrary large threshold as equivalent). This 
simple numerical method reveals the bifurcation scenario 
of the period-doubling route to chaos with the occurrence 
of several periodic windows in between, depicted by the 
red line. For comparison, we also show the average vari- 
ance of degrees of the corresponding networks. It is clear 
that this variance is in general larger in the chaotic re- 
gion. Note that when an incommensurate sampling fre- 
quency is chosen, the associated recurrence-based phase 
space network could become disconnected and be divided 
into several components rather than a regular single net- 



work at some values of parameter c in periodic region, 
inevitably leading to large variance in degree which ap- 
pears in even low-order periodic regions of the bifurcation 
diagram. Hence, we use a moving average with window 
width 3 for consecutive values to smooth the curve of the 
calculated values for variance of degrees in the blue line 
in Fig.[§b). 



Degree distribution 
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FIG. 6. (a) Color coded representation of vertex degree 
on the attractor for the chaotic Rossler system with c = 
18. (b) Zoom-in area in the dashed box of (a). 



B. Local network characteristics 

1. Local vertex degree 

Rather than simply looking at the global properties 
of the network, we also examine local vertex properties 
which we believe can provide more detailed information 
that is otherwise buried in the average geometry of the 
attractor. The first local property we study is the lo- 
cal vertex degree which is the number of neighbors of a 
given vertex v in the network. The degree of a vertex for 
the traditional recurrence networks which are based on 
an appropriate choice e reflects exactly the local recur- 
rence rate. When it comes to the phase space network, 
we should recall that the average degrees are 8 due to 
the specific construction methods we use to obtain the 
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network. Thus the local degrees depend closely on den- 
sity gradients of the phase space attractor. On the one 
hand, in the region where the density of the phase space 
is homogenous the corresponding structures of the net- 
work retain the homogeneity from the attractor and local 
degrees for the vertices will be around 8. That is, every 
point in these regions is connected to its 4 nearest neigh- 
bors and at the same time gains (on average) 4 extra links 
from points which take it as a nearest neighbor. One the 
other hand, if points are located in the regions of large 
gradient in the phase space density, some of these points 
may gain more than 4 extra links besides their 4 nomi- 
nal nearest neighbors, resulting in the lack of neighbors 
for some other points in the same region. These regions 
will have a heterogeneous degree distribution. Long pe- 
riodic signals follow a uniform distribution of low dimen- 
sion while the trajectories of chaotic signals tend to be 
trapped in the vicinity of the unstable periodic orbits 
(UPO) that lead to large gradients in the phase space 
density. One may refer to color-coded representation of 
the vertex degree on the attractor for the chaotic Rossler 
system with c = 18 in Figs. |6^a) and (b) for details. 
Overall, we get the heterogeneously distributed degrees 
in the attractor in Fig. |6](a). Typically, there are UPOs 
in the chaotic attractor. The trajectory will approach a 
UPO along its stable manifold, stay close to the UPO for 
several cycles and then be ejected away until it is cap- 
tured by the stable manifold of another UPO, leading to 
density variation. 
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FIG. 7. Local heterogeneity in terms of degree variance 
for the three orbits. 



Clustering coefficient distribution 
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To provide a more quantitative analysis, we introduce 
the local heterogeneity for a specific part of phase space 
which is defined by the variance of degrees of all points 
located within a ball of appropriate radius r in the phase 
space. By zooming in a particular region of the attractor 
showed in Fig. |6jb), we check three orbits representing 
trajectories in different situations with local heterogene- 
ity plotted in Fig. [7] Orbit 1 is in a low-density region 
while orbit 2 is located in the region close to a domi- 
nant UPO with a high density in this chaotic attractor. 
We observe a heterogeneous degree distribution in the 
nearby trajectory along orbits 1 and 2, because both of 
the two orbits are located in regions with higher gradient 
of attractor density. In comparison, we get a homogenous 
degree distribution in the region along orbit 3, which is 
located in the part of the attractor with lower gradient in 
density. If the radius r is chosen too large, the ball will 
cover a larger part of the attractor and local heterogene- 
ity will be removed. Hence we find that the three curves 
in Fig. [7] tend to converge to a fixed value as r increases. 
Moreover, we can address the diverse variance of degrees 
in Fig. [5] for different dynamical types. Since the attrac- 
tor becomes more heterogeneous as the dynamics changes 
from low-order periodic to high-order periodic and finally 
chaotic, the values for the degree variance will be larger. 




FIG. 8. (a) Color coded representation of the ver- 
tex clustering coefficient on the attractor for the chaotic 
Rossler system with c = 18. (b) Zoom-in area in the 
dashed box of (a). 



2. Clustering coefficient of a specific vertex 

The second property we study is the clustering coef- 
ficient of a specific vertex which quantifies the relative 
density of connections between its neighbors and thus 
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measures the network at a relatively large scale compared 
to the local degree. Figures [8]ja) and (b) show color rep- 
resentations of the vertex clustering coefficient on the 
attractor for the chaotic Rossler system with c = 18, 
from which we observe the heterogeneously distributed 
local clustering coefficients. One might expect a simi- 
lar effect for epsilon-recurrence networks within a certain 
range where a node in a dense region of the phase space 
will naturally have a larger clustering coefficient and vice 
versa. However, for recurrence-based phase space net- 
work, shown here in Figs. [8]ja) and (b), the values of 
the clustering coefficient for some nodes located in the 
denser regions of the attractor may be small. For exam- 
ple, clustering coefficients for the nodes located on the 
nearby orbit next to orbit 2 are small, while the values 
for nodes located in the boundary region of the attrac- 
tor can be large. In a homogeneously filled region, things 
become simpler and we can see points with relatively uni- 
form and higher clustering coefficients. However, when it 
comes to a region with large density fluctuations, we have 
to consider the interactions between several orbits which 
are close to each other. A strong connectivity between 
neighbors of the same point can be obtained only if its 
neighbors tend to lie in different orbits, because we have 
excluded links within the same orbit when constructing 
the phase space networks^. We can conclude that the 
local clustering coefficient contains the complexity of the 
interaction between orbits in the phase space and can 
provide complementary information on the system. 



3. Betweenness centrality of a specific vertex 



Betweenness centrality distribution 
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FIG. 9. (a) Color coded representation of the logarithm 
of vertex betweenness centrality, log(b„ + 1), on the at- 
tractor for the chaotic Rossler system with c = 18. (b) 
Zoom-in area in the dashed box of (a). 



The third property we study is the betweenness cen- 
trality of a specific vertex which quantifies the frequency 
with which shortest paths (between a randomly chosen 
pair of nodes) pass through the vertex. Thus it is more 
complicated because it measures the network at an even 
larger scale compared to the clustering coefficient and 
the local degree. What we are most interested in are the 
points with a large betweenness centrality which are of 
importance for many shortest paths on the network. As 
shown in Fig. [9] the nodes with high betweenness cen- 
trality may correspond to a region with low phase space 
density between such regions because of an increasing 
number of shortest path between such regions. 

In Fig. 10 we show the relationship between vari- 
ous local network properties. In Fig. flOta), the log- 



betweenness follows a bell-shaped distribution in which 
most of the points have around 10 4 — 10 5 shortest paths 

and (b) which shows 
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passing through. From Figs, 
the relationships of vertex betweenness centrality with 
local degree and local clustering coefficient, points with 
low betweenness (below 10 2 ) have relative low degree and 
high clustering coefficient. In both relationships, how- 
ever, we find a broad range of values of betweenness for 
most of points with respect to a fixed value of either de- 
gree or clustering coefficient. 
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Logarithm of vertex betweenness centrality, log(b v +1) 

FIG. 10. (a) Histogram for the logarithm of vertex 
betweenness centrality, log(&„ + 1). (b) Relationships 
between local degree, k and logarithm of betweenness 
centrality, log(6u + 1). (c) Relationships between local 
clustering coefficient, C v and logarithm of betweenness 
centrality, log(6„ + 1). 



As a consequence of the above observations, the com- 



plex network scheme can be used as a multi-scale tool to 
describe properties contained in time series from global 
to local scales. Moreover, different local vertex properties 
contain information on different aspects from the geom- 
etry of the attractor. 



III. EFFECT OF NOISE 

We now extend the multiscale characterization of com- 
plex networks obtained from time series to systems cor- 
rupted by noise. The first question here is: how will 
the noise affect the super-family phenomenon shown in 
Ref. [24t* To look for the answer, we systematically calcu- 
late the motif ranking for data from a different dynamical 
system corrupted by different levels of white Gaussian 
noise. The results are shown in Table |U When check- 
ing the first column in Table [I] which shows results for 
signals without noise, we obtain the superfamily phe- 
nomenon which is consistent with the results shown in 
Ref. [2H We keep the embedded parameters the same 
as for signals which are free of noise, because we want 
to investigate the effect of noise on the network struc- 
ture rather than that of the embedding parameters. Net- 
works from both continuous dynamics and map data ex- 
hibit superfamily phenomena that we can identify as pe- 
riodic, chaotic, hyperchaotic flow and maps by using the 
phase space method. For example, motif ranks for net- 
works from the periodic Rossler system are ABCFDE. 
Some networks completely lack motif D or E because 
the frequencies of the two are statistically insignificant 
for those networks. Chaotic flows such as Chua, Lorenz 
and chaotic Rossler systems with c = 18 exhibits the 
same motif rankings ABCDFE. The motif ranking for 
Mackey-Glass flow is ADBCEF which is similar to pe- 
riodic systems corrupted by noise because the system is 
in a high-dimensional state. For map data, the three 
chaotic maps (Henon map, Ikeda map and logistic map) 
all have motif ranks ABCFDE, while the hyperchaotic 
generalized Henon map and folded-tower map data have 
the same subgraph ranks ABCDFE and hence belong 
to one superfamily. 

We find that the superfamilies become less distinct 
when we gradually increase the amplitude of noise - 
once the noise becomes too strong, motif superfamilies 
lose their discriminatory power, as expected. By apply- 
ing an intermediate amount of noise we can see the switch 
in the sequence of motifs which shows the transitions 
from pure signals with no noise to signals with levels of 
noise. The rank of nontransitive motifs of type D rep- 
resenting a subgraph with 4 nodes and 3 edges which 
connect one node to the other three nodes is increasing, 
indicating that the system is higher-dimensional, or the 
distribution of nodes is more heterogeneous. By checking 
the last column in Table [I] we see results for signals with 
very significant noise (OdB noise). Here we observe that 
the networks are divided into two groups by motif ranks: 
all the flow data have motif ranks ADBCEF , while all 



the map data have motif ranks ABDCEF. Of course, for 
more moderate levels of noise, the power of the motif su- 
perfamily classification to correctly distinguish different 
dynamics improves. Finally, we observe the possibly sur- 
prising result that adding a higher level of noise actually 
helps to distinguish between maps and flows. Although 
this distinction is artificial (all the data we consider are 
by definition discretely sampled), it is effectively a mea- 
sure of the relative sampling rate of the dynamics and 
the discrete noise source. Systems with smooth dynam- 
ics and discontinuous dynamical systems become easy to 
distinguish when details of dynamical types are obscured 
by noise. 

We also examine the impact of noise on the standard 
measurements of the network including degree distribu- 
tion, average path length and clustering coefficient. For 
the global degree distribution that we have discussed 
above, we infer from the network construction methods 
that the average degree remains 8. In the presence of 
noise as shown in Fig. |11| however, the bell shape de- 
gree distribution loses its symmetry and becomes a skew 
distribution with mode 4. 

In Fig. [l2j we show how the average path length 
changes the network size for a particular realization in 
the presence of noise of different levels. The scaling be- 
havior for the average path length has been shown in the 
previous section. In comparison, for a particular run, 
only when the length of the signal is long enough can it 
contain sufficient phase space information, for example 
enough UPOs, to characterize chaos. Hence the aver- 
age path length is larger than the expected scaling value 
for networks from signals with length smaller than 7000. 
This does not affect the result, however, as we only try 
to explore the impact of noise. In the presence of a small 
amount of noise, the average path length decreases a lit- 
tle. When we increase the noise level, we still see the 
scaling curve increasing, but the average paths length 
becomes much smaller than the original one. 

Figure fl3] shows the clustering coefficient distribution 
for the same set of data with different noise levels. The 
phase space network exhibits a large clustering coefficient 
for the chaotic Rossler system in the absence of noise. We 
find the bell shape distribution for clustering coefficients 
with the maximum values around 0.5. Adding noise re- 
sults in the left to right shift in the distribution which 
means the noise tends to destroy the clustering struc- 
tures in the network. 

In the above discussion, we examined the impact of 
noise on the network properties from a global view. The 
existence of multiscale descriptions given by network the- 
ory will provide us with more detailed information by 
measuring the local vertex properties as well. In the next 
section we will demonstrate that this method can also be 
usefully applied to experimental time series data. In this 
case the experimental data comes from a "toy" physi- 
cal experiment. Nonetheless, the methods could also be 
applied to any system of particular interest. 
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TABLE I. Motif ranking for system corrupted by different levels of noise. 



System 


Dynamical type 


No noise 


30dB 


20dB 


lOdB 


OdB 


Rossler c — 6 


period = 2 


ABCFD 


ABDCEF 


ADBCEF 


ADBCEF 


ADBCEF 


Rnsslpr c — 1 9 


nprind — ^ 


ABCF 


ABDCFE 


ABDCEF 

ills I y V. 1 J I 


ADBCEF 


ADBCEF 


Rossler c — 8.5 


period— 4 


ABCFD 


ABDCFE 


ADBCEF 


ADBCEF 


ADBCEF 


Rossler c = 12.6 


period= 6 


ABCDFE 


ABDCFE 


ABDCEF 


ADBCEF 


ADBCEF 


Rossler c = 8.7 


period= 8 


ABCDFE 


ABDCFE 


ADBCEF 


ADBCEF 


ADBCEF 


Rnsslpr p — Q 


ph a nt i p fl ow 


ABCFDE 


ABDCFE 

a \ l j I y v_ i i j 


ADBCEF 


ADBCEF 


ADBCEF 


Rossler c — 18 


chaotic flow 


ABCDFE 


ABCDFE 


ABDCEF 


ADBCEF 


ADBCEF 


Lorenz 


chaotic flow 


ABCDFE 


ABCDFE 


ABDCFE 


ADBCEF 


ADBCEF 


Chua 


chaotic flow 


ABCDFE 


ABDCFE 


ADBCEF 


ADBCEF 


ADBCEF 


Mackey-Glass 


chaotic flow 


ADBCEF 


ADBCEF 


ADBCEF 


ADBCEF 


ADBCEF 


Logistic 


chaotic map 


ABCFDE 


ABDCFE 


ABDCEF 


ABDCEF 


ABDCEF 


Henon 


chaotic map 


ABCFDE 


ABDCFE 


ABDCEF 


ABDCEF 


ABDCEF 


Ikeda 


chaotic map 


ABCFDE 


ABDCEF 


ABDCEF 


ABDCEF 


ABDCEF 


hyperHenon 


hyperchaotic map 


ABCDFE 


ABCDFE 


ABDCFE 


ABDCEF 


ABDCEF 


hyperTower 


hyperchaotic map 


ABCDFE 


ABDCFE 


ABDCFE 


ABDCEF 


ABDCEF 



IV. ANALYSIS OF TIME SERIES FROM CLARINET 
DATA 

As a real-world example, we apply the network method 
to experimental data of sustained tones voiced on a stan- 
dard Bb clarinet over most of the dynamic range of the 
instrument — from E 3 to F 6 in standard scientific nota- 
tion. 

A similar application was initially described in Ref. 1261 
The clarinet tones are produced by self-sustained oscil- 
lations coupling with nonlinear interactions between the 
flow and the mechanical response of the reed. By mak- 
ing use of the same complex network method for several 
sets of time series data from the clarinet, the preliminary 
result which has been reached is that as the data belongs 
to superfamilies which are most similar to chaotic or hy- 
perchaotic dynamics, the pure sustained clarinet tones 
could be characteristically aperiodic, and consistent with 
chaotic or hyperchaotic dynamics. However, the dynam- 
ics in the clarinet tones which we observed could be af- 
fected by either the performer's skills or the quality of 
the clarinet. 

In this paper, to extend the same analysis method on 
clarinet tones, we make use of another set of clarinet data 
from the stationary period of 23 distinct notes individu- 
ally recorded at a sampling frequency of 44.1kHz. After 
being smoothed by an appropriate window width, each 
data was down-sampled to a level with approximately 25 
samples per cycle so that it covers more or less the same 
number of orbits in the phase space. The time delay was 
chosen to be the first minimum of mutual information 
(typically between 2 and 8) and the embedding dimension 
was 10 which we verified to be significantly larger than 
the estimated correlation dimension of all data (between 
1 and 3) . By linking each vertex with its 4 nearest neigh- 
bors, the phase space networks have been constructed 
from the embedded time series. 



A single motif ranking, ABDCEF dominates. In all 
but two cases, ABCDFE or ABDCEF is the appropri- 
ate motif ranking. The two exceptions are Cg and G3. 
The note C6 corresponds to one at the very top of the 
range of the second register (the so-called clarion regis- 
ter) and is obtained with the very shortest possible air 
column length, and the fundamental harmonic removed. 
This note therefore represents an extremity of the range 
of the instrument. Of course, higher notes are possible 
(into the altissimo register — such notes are obtained 
by removing still further harmonics), however, these are 
produced with much less regular fingering patterns and 
require tighter control of the reed by the musician - 
to ensure that those lower harmonics are correctly re- 
moved. As far as we know, the note G3 is not peculiar 
— it is one of the lowest notes of the bottom (chalumeau 
register), but not the lowest. The predominant motif 
ranking which we observe corresponds to that for chaotic 
and periodic flows with moderate noise. The large path 
length and diameter are consistent with periodic flows 
exhibiting small to moderate noise. The moderate clus- 
tering coefficient and large bctweenncss arc similarly con- 
sistent with a large degree of spatial heterogeneity in the 
underlying attractor. In fact, more local information is 
contained in the local scale rather than the average be- 
tweenness properties. 



While we have not been able to provide a definitive 
answer to the question of whether clarinet tones are 
"chaotic" , the evidence we have obtained is certainly not 
inconsistent with that proposition. Of course, the issue of 
whether the intonation could be generated by a noisy pe- 
riodic orbit remains unaddressed. To resolve this we will 
need to resort to surrogate data methods ^Ht-HD Nonethe- 
less, it is encouraging that consistent results have been 
obtained across the full spectrum of experimental data. 
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TABLE II. Motif ranking for experimental data of sustained tones voiced on a standard B\> clarinet 



Notes 


Motif ranking 


Mean be 


mean cc 


average path 


diameter 


1 


ARDCEF 


4 07068fip4-oo5 


fl 3Q31 5 


41.7109 


86 


Fq 


ARDCEF 


4 fi47970p+005 


o 357ifi 


47 4843 


96 




ADRCEF 


3 598no8p+oo5 


31 33Q 


3fi 9837 


75 


A, 


ARDCEF 


7 344845p4-005 


fl 3Q7Q4 


74.0342 


154 


Rq 


ARDCEF 


i oofifvwp+oos 


0.37172 


71.5472 


144 


C,A 


ARDCEF 


1 337349p+01 9 


fl 4fi934 


83 5449 


169 


U4 


ARDCEF 

f\-±J 1/ k_' J-J J. 


5 555790P+005 


fl 40947 


5fi 5fi1 8 


113 




ARDCEF 


4 9735fifip-l-01 5 


fl 39809 


98.4771 


200 


F/ 
£ 4 


ARDCFE 


4 994994p+008 


fl 41 fi9fi 


fil 991 3 


125 




ARDCEF 


4 8fifi59fip4-oo7 


fl 37SR7 


74.1277 


150 


A* 


ARDCEF 


fi fi9i fi53p+oo5 


3745fi 


fi5 0096 


131 


±J>4 


ARDCEF 


fi 973fi7fip+01 fi 


fl 41 S9Q 


95 0875 


192 


^5 


ARDCEF 


fi 4Q1 895p+005 


fl 3233Q 


fi5 0739 


132 


J-A5 


ARDCEF 


5 083065p-i-oi q 


0.40712 


1 1 9 3084 


240 




ARDCEF 


K 77QS99P+005 


38^8(1 

U .OOJOU 


58 7799 


119 


-L 5 


ARDCFE 


7 343508p+oi 2 


0.44612 


79.1254 


160 


G 5 


ABDCFE 


1.531236e+007 


0.4283 


57.275 


114 


A 5 


ABDCFE 


2.127506e+007 


0.4472 


67.9874 


137 


B 5 


ABDCEF 


1.074712e+006 


0.3806 


86.8553 


175 


c 6 


ABCDFE 


1.394829e+030 


0.56559 


100.6292 


205 


D 6 


ABDCEF 


2.350170e+020 


0.41137 


94.4897 


191 


E 6 


ABDCFE 


1.877442e+019 


0.50417 


107.3988 


217 


F 6 


ABDCEF 


1.663571e+013 


0.31788 


114.8136 


232 



V. CONCLUSIONS 

The purpose of this paper has been to provide a more 
thorough investigation of the characteristic properties of 
complex networks obtained from times series, specifically 
using the phase space method of Xu et aipS. To do so 
we have not only provided a more detailed exploration of 
the effect of noise contamination on the motif superfamily 
methods, but we have also extended the analysis of cor- 
responding networks to include micro- and meso- scale 
properties as well as the usual global summary statis- 
tics. This approach has provided new insights as we have 
been able to connect dynamical properties of the under- 
lying system to variations in complex network properties 
of individual nodes. In particular, average path length, 
global clustering coefficient and degree distribution de- 
fined in the global network scale can be used to distin- 
guish systems of different dynamical types. Local degree, 
clustering coefficients and betweenness centrality contain 
information on the spatial fillings of the phase space and 
are close to invariant objects such as UPOs. The specific 
vertex properties can provide more detailed information 
about the local attractor geometry in the phase space. 

We demonstrated that the method can also be usefully 
applied to real experimental data without being over- 
whelmed by the effect of noise. In this case, the natural, 
but still rather obvious conclusion is that clarinet sound 
production is consistent with either pseudoperiodic chaos 
or a noise contaminated periodic orbit. Nonetheless, this 



is a useful conclusion as methods exist to distinguish be- 
tween these two cases, and the application of such meth- 
ods with the measures described here as test statistics 
will be straightforward 2211 ^. 
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